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Abstract 

The properties of the partition function zeros in the complex temperature plane (Fisher zeros) 
and in the complex Q plane (Potts zeros) are investigated for the Q-state Potts model in an arbitrary 
nonzero external magnetic field H q , using the exact partition function of the one-dimensional model. 
The Fisher zeros of the Potts model in an external magnetic field are discussed for any real value 
of Q > 0. The Potts zeros in the complex Q plane for nonzero magnetic field are studied for the 
ferromagnetic and antiferromagnetic Potts models. Some circle theorems exist for these zeros. All 
Fisher zeros lie on a circle for Q > 1 and H q > except Q = 2 (Ising model) whose zeros lie on the 
imaginary axis. All Fisher zeros also lie on a circle for any value of Q when H q = (except Q = 0, 
1 and 2) or H q = — oo (except Q = 1, 2 and 3). All Potts zeros of the ferromagnetic model lie on 
a circle for H q > 0. When H q = or — oo, all Potts zeros lie on a circle for both the ferromagnetic 
and antiferromagnetic models. All Potts zeros of the antiferromagnetic model with H q < also lie 
on a circle for (x + l) -1 < a < 1, where a = e@ J and x = e /3Hq . It is found that a part of the Fisher 
zeros or the Potts zeros lie on a circle for the specific ranges of H q . It is shown that some Fisher 
or Potts zeros can cut the positive real axis. Furthermore, the Fisher zeros or the Potts zeros lie 
on the positive real axis for the specific ranges of H q . The densities of zeros are also calculated 
and discussed. The density of zeros at the Fisher edge singularity diverges, and the edge critical 
exponents at the singularity satisfy the scaling law. There exists the Potts edge singularity in the 
complex Q plane which is similar to the Fisher edge singularity in the complex temperature plane. 

PACS numbers: 05.50.+q; 64.60.Cn; 75.10.Hk; 02.10.Ox 
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I. INTRODUCTION 

The Q-state Potts model □ DDDD is a generalization of the Ising (Q = 2) model. 
The Q-state Potts model exhibits a rich variety of critical behavior and is very fertile ground 
for the analytical and numerical investigations of first- and second-order phase transitions. 
The Potts model is also related to other outstanding problems in physics and mathematics. 
Fortuin and Kasteleyn 0, 3} have shown that the Q-state Potts model in the limit Q — > 1 
defines the problem of bond percolation. They |7[ also showed that the problem of resistor 
network is related to a Q = limit of the partition function of the Potts model. In addition, 
the zero-state Potts model describes the statistics of treelike percolation & and is equivalent 
to the undirected Abelian sandpile model The Q — \ state Potts model has a connection 



to a dilute spin glass 



jlpl |. The Q-state Potts model with < Q < 1 describes transitions in 



the gelation and vulcanization processes of branched polymers 

By introducing the concept of the zeros of the partition function in the complex magnetic- 
field plane (Yang-Lee zeros), Yang and Lee 3] proposed a mechanism for the occurrence of 
phase transitions in the thermodynamic limit and yielded a new insight into the unsolved 
problem of the ferromagnetic Ising model in an arbitrary nonzero external magnetic field. 



It has been known 
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17, 



18, 



19 . 120L l2l| that the distribution of the zeros of 



a model determines its critical behavior. Lee and Yang [13J also formulated the celebrated 
circle theorem which states that the partition function zeros of the Ising ferromagnet lie 
on the unit circle in the complex magnetic-field (x = e 1311 ) plane. However, for the ferro- 
magnetic Q-state Potts model with Q > 2 the Yang-Lee zeros in the complex x plane lie 
close to, but not on, the unit circle with the two exceptions of the critical point x = 1 
{H = 0) itself and the zeros in the limit T = Q fi H Q Q. It has been shown 
24 , |3| that the distributions of the ferromagnetic Yang-Lee zeros for Q > 1 have similar 



properties independent of dimension. Recently, the exact results on the Yang-Lee zeros of 
the ferromagnetic Pott^ model have been found using the one-dimensional model 24 , 2|| . 



Mittag and Stephen W stndied thc^ng-Lcc Z eros of the three-state Potts ferromagnet in 
one dimension. Glumac and Uzelac |27J] found the eigenvalues of the transfer matrix of the 
one-dimensional Potts model for general Q. In particular, they have shown that the ferro- 
magnetic Yang-Lee zeros can lie on the positive real axis for Q < 1. Yang also calculated 
the Yang-Lee zeros of the one- dimensional Ising antiferromagnet, and showed that they lie 
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on the negative real axis 



28|. 



The first zero, which we call the edge zero, and its complex conjugate of a circular 
distribution of the Yang-Lee zeros of the Potts model cut the positive real axis at the 
physical critical point x c = 1 for T < T c in the thermodynamic limit. However, for T > T c 
the edge zero does not cut the positive real axis in the thermodynamic limit, that is, there 
is a gap in the distribution of zeros around the positive real axis. Within th is g ap, the 
free energy is analytic and there is no phase transition. Kortman and Griffiths |29| carried 
out the first systematic investigation of the magnetization at the edge zero, based on the 
high-field, high-temperature series expansion for the Ising model on the square lattice and 
the diamond lattice. They found that above T c the magnetization at the edge zero diverges 
for the square lattice and is singular for the diamond lattice. For T > T c we rename the edge 
zero as the Yang-Lee edge singularity. The divergence of the magnetization at the Yang-Lee 
edge singularity means the divergence of the density of zeros, which does not occur at a 
physical critical point. Fisher (^J proposed the idea that the Yang-Lee edge singularity can 
be thought of as a new second-order phase transition with associated critical exponents and 
the Yang-Lee edge singularity can be considered as a conventional critical point. The critical 
point of the Yang-Lee edge singularity is associated with a 3 theory, different from the usual 
critical point associated with the 4 theory. The crossover dimension of the Yang-Lee edge 
singularity is d c = 6. 

Fisher jlijj emphasized that the partition function zeros in the complex temperature plane 
(Fisher zeros) are also very useful in understanding phase transitions, and showed that for 
the square lattice Ising model in the absence of an external magnetic field the Fisher zeros in 
the complex y = e~^ J plane lie on two circles (the ferromagnetic circle y-pu = —l + \/2e %e and 
the antiferromagnetic circle Vafm = 1 + V2e t9 ) in the thermodynamic limit. In particular, 
using the Fisher zeros both the ferromagnetic phase and the antiferromagnetic phase can be 
considered at the same time. The critical behavior of the square-lattice Potts model in both 
the ferromagnetic and antiferromagnetic phases has been studied using the distribution of the 
Fisher zeros, and the Baxter conjecture for the antiferromagnetic critical temperature 



has been verified \3l 



341 ] . Recently, the Fisher zeros of the Q-state Potts model on square 



lattices have been studied extensively for integer 
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Q > 2 



Q, 0,0, EL 



numerical studies have shown 
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Exact 



47| that for self-dual boundary 
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conditions the Fisher zeros of the Q > 1 Potts models on a finite square lattice are located 
on the unit circle in the complex p plane for Re(p) > 0, where p = (e^ J — 1)/ \[Q. It has been 
analytically shown that all the Fisher zeros of the infinite-state Potts model lie on the unit 
circle for any size of square lattice with self-dual boundary conditions (Q, and the Fisher 
zeros near the ferromagnetic critical point of the Q > 4 Potts models on the square lattice 
lie on the unit circle in the thermodynamic limit Chen et al. j^J conjectured that 

when Q reaches a certain critical value Q C (L), all Fisher zeros for L x L square lattices with 
self-dual boundary conditions are located at the unit circle \p\ = 1. Later, Kim and Creswick 
|s4 | verified this conjecture and found that Q C (L) approaches infinity in the thermodynamic 
limit. The thermal exponent y% of the square- lattice Potts antiferro magnet has also been 



studied using the Fisher zeros near the antiferromagnetic critical point 

The Yang-Lee zeros in the complex x plane for real y have been extensively studied and 
well understood. However, the Fisher zeros in the complex y plane for real x are much 
less well understood than the Yang-Lee zeros in the complex x plane. The Fisher zeros in 

n 

an external magnetic field have been considered by Itzykson et al. |14( for the first time. 
They studied the movement of the Fisher zero closest to the positive real axis for the Ising 
model as the strength of a magnetic field (x > 1) changes. For nonzero magnetic field 
there is a gap in the distribution of the Fisher zeros of the Ising model around the positive 
real axis even in the thermodynamic limit, which means that there is no phase transition. 
Matveev and Shrock 0, Q studied the Fisher zeros of the two-dimensional Ising model 
in an external magnetic field {x > 1) using the high-field, low-temperature series expansion 
and the partition functions of finite-size systems. They observed that the thermodynamic 
functions and the density of zeros diverge at the edge zero when x > 1, and that the values 
of the critical exponents associated with the edge zero for x > 1 are nearly independent 
of x and satisfy the Rushbrooke scaling law approximately. We call the edge zero in the 
complex temperature plane for x > 1 as the Fisher edge singularity. Kim and Creswick 
extended the study of the Fisher edge singularity to the Q-state Potts model (Q > 3) on the 
square lattice. They found that the values of the edge critical exponents at the Fisher edge 
singularities are nearly independent of Q (Q > 2). The critical behavior at the Fisher edge 
singularity has also been studied for the Ising model and the three-state Potts model on the 
simple-cubic lattice using the high-field, low-temperature series expansions (^j). However, 
the Fisher zeros of the Q-state Potts model in an external magnetic field have never been 
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investigated for general Q except few integer values of Q. 

The partition function of the Q-state Potts model in the absence of an external magnetic 
field is also known as the Tutte dichromatic polynomial [3| or the Whitney rank function [^j] 
in graph theory and combinatorics of mathematics. One of the most interesting properties 
of the antiferromagnetic Potts model is that for Q > 2 the ground-state is highly degenerate 
and the ground-state entropy is nonzero. The ground states of the antiferromagnetic Potts 
model are equivalent to the chromatic polynomials [56] in mathematics, which play a central 
role in the famous four-color problem [57(. The partition function zeros in the complex Q 
plane (Potts zeros) of the Q-state Potts model have been studied both in mathematics and 
in physics. The Potts zeros at j3J = — 00, corresponding to T = for the antiferromagnetic 
model, have been investigated extensively to understand the chromatic polynomials and the 



ground states of the antiferro mag netic Potts model 
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^ Recently, the 
L x L square 



Potts zeros at finite temperatures have been studied for ladder graphs 

lattices , and the model with long-range interactions 0| . However, the properties of the 
Potts zeros for nonzero magnetic field have never been known until now. 

In this paper, we investigate the exact results on the partition function zeros of the 
ferromagnetic and antiferromagnetc Q-state Potts models in one dimension, and we unveil 
the unknown properties of the Fisher zeros and the Potts zeros in an external magnetic field. 
In the next section we obtain two master equations to determine the partition function zeros 
of the one- dimensional Potts model. In Sec. Ill we calculate and discuss the Fisher zeros of 
the Q-state Potts model for any value of Q and in an arbitrary magnetic field. In Sec. IV 
we study the Potts zeros in the complex Q plane of the ferromagnetic and antiferromagnetic 
Potts models in an arbitrary external magnetic field. 



II. THE Q-STATE POTTS MODEL 

The Q-state Potts model in an external magnetic field H q on a lattice G with N s sites 
and Nf, bonds is defined by the Hamiltonian 

Hq = -Jj2 s ( a i> a i) ~ H vl2 s ( a k,q), (l) 

(id) * 
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where J is the coupling constant (ferromagnetic model for J > and antiferromagnetic 
model for J < 0), indicates a sum over nearest-neighbor pairs, Oi = 1, 2, Q, and q is 
a fixed integer between 1 and Q. The partition function of the model is 

Z Q = E e"^, (2) 

where {o" n } denotes a sum over possible spin configurations and (3 = (ksT)" 1 . The 
partition function can be written as 

N b N s 

Z(a, x, Q) = E E V Q (E, M)a E x M , (3) 

where a = y~ l = e@ J , x = e^ Hq , i£ and M are positive integers < E < and < M < N s , 

respectively, and Qq(E,M) is the number of states with fixed E and fixed M. The states 

with E = (E = Nj,) correspond to the antiferromagnetic (ferromagnetic) ground states. 

The parameter Q enters the Potts model as an integer. However, the study of the Q-state 

Potts model has been extended to continuous Q due to the Fortuin-Kasteleyn representation 

of the partition function 0, 3] and its extension 

For the one- dimensional Potts model in an external field the eigenvalues of the transfer 

I — I ° 

matrix were found by Glumac and Uzelac 27]. The eigenvalues are X± = (A±iB) /2, where 



A = a(l + x)+Q-2 and B = -i^[a{l - x) + Q - 2] 2 + 4(Q - l)x, and A = a- 1 which is 
(Q— 2)-fold degenerate. The partition function of the one-dimensional model (N = N s = Nb) 
is given by 

Z N (a,x 1 Q) = \^ + \^ + (Q-2)\£. (4) 
When A + and A_ are two dominant eigenvalues, we have 

Z N ~ A? + X N (5) 



for large N. If we define A = 2C cosip and B = 2C sin?/', where C = J (a — l)(a + Q — l)x, 
then A± = C exp(±iip), and the partition function is 



l N = 2C N cos Nip. (6) 



The zeros of the partition function are then given by 



2k 4- 1 

V> = ^fc=-2^-T, fc = 0,l,2,...,iV-l. (7) 



In the thermodynamic limit the locus of the partition function zeros is determined by the 
solution of 

A = 2Ccos^, (8) 

where < ip < ir. In the special case Q = 2 the contribution by the eigenvalue Ao disappears 
from the partition function, Eq. (4), and the equation (8) determines all the locus even for 
finite systems. 

On the other hand, when A + and A are two dominant eigenvalues, we have 

Z N ~\» + (Q-2))$ (9) 
for large N. The partition function zeros are then determined by 

^± = (2-Q) 1 ^ eX p(#), (10) 

where 

= fc = ^, fc = 0,l,2,...,7V-l. (11) 

In the thermodynamic limit the locus of the partition function zeros is determined by the 
solution of 

aV + (Q- l)x - axA + (a - l)Ae^ - (a - 1)V^ = 0, (12) 

where < < 2ir. The equation (12) also determines the locus of the zeros when A_ and 
A are two dominant eigenvalues. 



III. FISHER ZEROS 

From the equation (8) the locus of the Fisher zeros is obtained to be 

(Q-2)0rcos2^-l)±22v/7 
y^>~ (Q-2)2 + 4(Q-l):rcos 2 ^ ' 1 ' 

where / = x cos 2 ip[(Q + x — l)(Qx — x + 1) — Q 2 x cos 2 if)]. The edge zeros of yi(ip) are given 

by 

(Q - 2)(x - 1) ± i2\x - 1\J(Q - l)x 

In the absence of an external field, H g — (x — 1), these edge zeros cut the real axis at the 
origin, 

y± = o, (15) 
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which corresponds to the T = ferromagnetic transition point. If / < 0, the zeros of yiijp) 
lie on the real axis. However, if / > 0, the zeros of yi(tp) lie on a circle 



JJi 



{i>) = y c + De ±w ^ (16) 



(17) 



in the complex y plane, where y c (the center of the circle) and D are given by 

_ 1 (Q-2) 2 (x-l) 2 -E(x+l) 2 
Vc ~ 2(Q - 2) (Q- 2) 2 (x - 1) + E(x + 1) 

and 

1 (Q-2) 2 (3x+l)(:r-l) + E(x + l) 2 
2(Q-2) (Q-2)*(x-l) + E(x + l) ' 1 ' 

E is defined by 

E = (Q - 2) 2 + 4(Q - l)x, (19) 

the argument 6 is given by 



cos^(^) = 



(Q - 2) (x cos 2t/> - 1) 



(20) 



.(Q - 2) 2 + 4(Q - l)xcos 2 V " 
and the radius r of the circle is 

r=\D\. (21) 



The one point of the circle yi(V'), 



^ = I or9 = = (22) 

always lies on the real axis. The point yi(f) lie on the positive real axis for Q < 2 and on 
the negative real axis for Q > 2. 

On the other hand, the equation (12) gives the second locus of the Fisher zeros 

_ e ^-(x + l)e^ + x 

Vm ~ e ** + (Q-2)e»-{Q-l)x- (23) 

The two points of y2{4>), 

2/2 (0) = (24) 

and 

can be lie on the real axis. From |A±| = |Ao| we also obtain 

y- = {Q -~iy I + i < 26 > 
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A. Q > 1 



In the special case Q = 2 all the Fisher zeros lie on the imaginary axis, and they meet the 
real axis at the origin only for x — 1. At x — 1 (H q = 0), the loci yi(il>) and 2/2(0) become 
the identical locus circle |4a| 

^> = -oVoV < 27 > 

for any value of Q except Q = 0, 1, and 2. On this circle three eigenvalues have the same 
magnitude 

IX |2 IX |2 IX |2 (Q-1)(Q-1-2CQS P ) + 1 

2(1 — cos p) 

This circle cuts the real axis at two points 

y(p = 0) = (29) 

and 

The point y(p = 0) is the T = ferromagnetic transition point. The point y(p = tt) has no 
physical meaning for Q > 2, but it may correspond to an antiferromagnetic transition point 
for Q < 2 because the physical interval is 

1<2/<oo (oo>T>0) (31) 

for antiferromagnetic interaction J < 0. For x > 1 > 0) all the Fisher zeros lie on the 
circle yiijp) (f > for Q > 1 and x > 1) , whereas for x < 1 (if g < 0) they are located on 
the loop y 2 (4>)- 

Figure 1 shows the locus of the Fisher zeros in the complex y plane of the three-state 
Potts model for x = 1, 2 and 3. For Q = 3 the center y c and the radius r of the circle yi{ip) 
are given by 

y c =-l, r = l (x = l), (32) 
19 20, 1 

vc = -w> r= u (x = 2) ' (33) 

and 

33 35 , _ . , „ , 
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For x = 1 two edge zeros y± cut the real axis at the origin which is the T = ferromagnetic 
transition point. However, as x increases, they move away from the origin, and there is 
a gap in the distribution of zeros, centered at 9 = 0, that is, the density of zeros is zero, 
g{9) = 0, for \9\ < 9 . The edge angle 9 (= 9(ip = 0)) is given by 



cos 6» = 



(Q-2)(s-l) 
E Vc . 



(35) 



and the edge zeros are expressed as 

y ± =y c + De ±i6 °. (36) 
For example, for the three-state Potts model, the edge zeros are located at 

y ± = (x = l), (37) 
I/ ± = l(l±4i) (x = 2), (38) 

and 

y± = ^-(2±4VQt) (x = 3). (39) 
On the circle yi(ip), the density of zeros g(9(ip)) is given by 

I sin 6* I , , , 

9{0) = 1 1 , (40) 

27r v /(g - 2)(x - 1) - E(y c + D cos 0) 

where V is defined by 

v = D[(Q-2y + 2(Q-l)(Q-2)(x+l)] 

[(Q - 2) - 2(Q - l)(y c + D cos 9)]yJ(Q - 2)(x + 1) + (Q - 2) 2 (y c + Dcos9) ' 

For x = 1 the density of zeros reduces to a simple form 

a(9) = 1 ^~ 2 I (AO) 

y{ ' 27r(Q-2) 2 + 2(Q-l)(l-cosfl)' 1 J 



which gives 



9 <°> - s vh\ < 43 » 



at the edge zeros y± = 0, showing a first-order transition. 

On the other hand, for x > 1 the density of zeros near the edge zeros is expressed as 

9(0 ~ ^o) ~ 1 (44) 

y/y(e)- y± (e Q )' 
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that is, the density of zeros diverges at the Fisher edge zeros y± for x > 1. In this case, 
the Fisher edge zero is called the Fisher edge singularity because of the divergence of the 
density of zeros. The edge critical exponents a e , (3 e and j e associated with the Fisher edge 
singularity are defined in the usual way, 

C e ~(y-y±r a % (45) 

m e ~(y-y±f% (46) 

and 

Xe~(y-y±r je , (47) 

where C e , m e , and Xe are the singular parts of the specific heat, magnetization, and sus- 
ceptibility, respectively. The density of zeros near the Fisher edge singularity is also given 
by 

g(9~9 )~(y-y ± ) 1 - a e. (48) 

The Fisher edge singularity is characterized by the edge critical exponents a e = |, [3 e = — |, 
and 7 e = |. These values satisfy the Rushbrooke scaling law \2oL l47j 

a e + 2f3 e + 7e = 2. (49) 

Figure 2 shows the locus of the Fisher zeros in the complex y plane of the Q — | state 
Potts model for x = 1 and x = 3. For Q = | the values of the center y c , the radius r, and 
the edge zeros y± of the circle yi(ip) are 

y c = 2, r = 2, y± = (x = 1) (50) 

and 

2/c = ^, r = ^, y ± = l(-l±2V6i) (x = 3). (51) 

The overall behavior in the figure 2 is similar to that in the figure 1. However, there is a big 
difference in the location of the point yi(ip = | or 6 = ir), Eq. (22), of the circle ?/i(6 l ('0)). 
The point yi(7r) lies on the negative real axis for Q > 2, whereas it lies on the negative real 
axis for 1 < Q < 2. For 1 < Q < 2 and x > 1 we have 

2/i(tt) > 2, (52) 

and the point yi(7i) may correspond to an antiferromagnetic transition point. 
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In the limit H q — > oo (x — > oo) the positive field H q favors the state q for every site 
and the Q-state Potts model is transformed into the one-state model. The Fisher zeros are 
determined by 

Z(y, x - oo, Q) ~ n Q( E ^ M = N *)V~ E = (53) 

E=0 

for any Q. Because Qq(E, M = N s ) = 1 for E = Nb and otherwise, Eq. (53) is 

y- Nb = 0. (54) 



As x increases, \y\ for all the zeros increases without bound 

Now we turn to the distributions of Fisher zeros for x < 1 (H q < 0) which are completely 
determined by the loop 2/2(0)- The locus cuts the real axis at two points 2/2(0) and 2/2 (vr). 
The point 

2/2(0) = (55) 

is the T = ferromagnetic transition point. For Q > 3 the point 2/2(71") lies on the negative 
real axis, and it has no physical meaning. Figure 3 shows the locus of the Fisher zeros in 
the complex y plane of the three-state Potts model for x = \, which gives 

2/ 2 (tt) = -3. (56) 

For 1 < Q < 2, the point 2/2 (tt) lies on the positive real axis, and it may be considered an 
antiferromagnetic transition point because 

2/ 2 (tt) > 1. (57) 

Figure 4 shows the locus of the Fisher zeros in the complex y plane of the Q — | state Potts 
model for x = \. For Q = | and x = \ we obtain 

1 2 

V2^) = -. (58) 

For 2 < Q < 3, the value of 2/2 (tt) diverges at 

a* = (59) 

For 2 < Q < 3 and xi < x < 1, the point 2/2 (tt) lies on the negative real axis, and the 
distributions of the Fisher zeros show similar behavior to that in Fig. 3. However, for 
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2 < Q < 3 and x < x, 1/2(71") > 1, and the distributions of the zeros are similar to that in 
Fig. 4. 

In the limit H q — > — 00 (x — > 0) the partition function of the Q-state Potts model, Eq. (4), 
becomes 

Z(a, x = 0, Q) = (a + Q - 2) N + (Q - 2) (a - 1)*. (60) 

For an external field H q < 0, one of the Potts states is supressed relative to the others, and 
the symmetry of the Hamiltonian is that of the (Q — l)-state Potts model in zero external 
field. Therefore, the partition function Z(a,x — 0, Q + 1) of the (Q + 1) state Potts model 
in the limit H q — ► — 00 is the same as the partition function 

Z(a, x = 1, Q) = (a + Q - l) N + (Q - l)(a - 1)* (61) 

of the Q-state Potts model in the absence of an external magnetic field. As x decreases from 
1 to 0, the Q- state Potts model is transformed into the (Q — l)-state Potts model in zero 

n 

external field j42(. At x = 0, the Fisher zeros of the Q-state Potts model lie on a circle 

for any value of Q except Q = 1, 2, and 3. The zeros lie on the imaginary axis for the 
three-state Potts model at x — 0. 

B. Q < 1 and x > 1 

For a; < £2, where 

x 2 = ^ ( 63 ) 
the locus consists of the loop 2/2(0) (0* < < 27r — 0*) and the line yi(i[)) (0 < -0 < tp*) 
inside the loop 1/2- The line yi(ip) lies on the real axis between = y\(ip*) and ?/_ = ?/i(0) 
(y* < 2/- < 0), and meets with the loop 1/2(0) at the point 

V* = 2/2(0*) = 2/2(2vr - 0*), (64) 

where |A + | = |A_| = |A |. The loop 2/2(0) also cuts the real axis at the point 2/2 (tt) (> !)• 
Figure 5 shows the locus of the Fisher zeros for Q = \ and x = | which has 

2 20 
y„ = -2, y_ = 1 - -j= = -0.155, 2/2W = — . (65) 
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At x = X2, y* = — oo, and the line yi(tp) lies on the real axis between — oo and 2/- (< 0). 
In the region xi < x < £3, where 

13 = r75- (66) 

the line yi(ip) lies outside the loop 1/2(0)- The line y\ lies on the real axis between —00 and 
2/_ (< 0) and between 2/* (> 1) and 00. The loop 1/2(0) cuts the real axis at the point 2/2(71"), 
and meets with the line yi(tp) at 2/* (1 < 2/2(71") < y*). Figure 6 shows the locus of the Fisher 
zeros for Q = \ and re = 3, from which we obtain 

V- = ^(3 - 2>/6) = -0.506, y 2 (7r) = 2, y* = 4. (67) 

At x — x 3 , 0* = n, the other edge zero y + appears, and the loop 2/2(0) shrinks to the 
point 



y + = y, = = Q + (Q 2 _^. (68) 

For x > X3, the only locus is the line y\(ip) on the real axis between —00 and 2/- (< 0) and 
between y + (> 1) and 00. At Q — 0, two edge zeros are simply expressed as 

y± = — 2 — ^ ^ 

for x > 1, and all the Fisher zeros still lie on the real axis. 



C. Q < 1 and x < 1 

For x > 1 — Q, the locus consists of the loop 2/2(0) (—0* < < 0*), the line 2/1 (VO 
(0 < -0 < -00 and 7r — ip < -0 < -0*) on the real axis between 2/* = yi(ip*) and 2/- = 2/1 (0) 
(0 < 2/* < 2/- < 1)? an d the circle 2/1 ("0) ("00 < "0 < tt — -0 O )- The loop 2/2 meets with the line 
2/1 at the point 

y* = 2/2(0*) = 2/2(-0*)> (70) 

where |A + | = |A_| = |A |. The circle 2/i("0) cuts the real axis at two points y\ and 2/i(f) 
(2/* < 2/i < 2/- < 2/i (f ))> where the points 2/? and 2/1 (f) are given by 

/ / \ < 1 \ ^ n 2(Q-2)x(x-l) 

vi = 2/1(^0) = 2/1 (tt - rh) = y c + d = {Q _ mx _ l) + E{x + l) (7i) 

and 
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respectively. Similarly, the loop also cuts the real axis at two points 1/2(0) and (0 = 
2/2(0) < y*). For example, figure 7 shows the locus for Q = | and x = ^. For these values 
we obtain 

». = 4rf=^»- = ^(» + 6 ) /D= 0.947, tt (I)=H, (73) 
and the center y c and the radius r for the circle 

Vc = 165' ' = 165- (74) 

At rr = 1 — Q, the circle yi(ip) shrinks to the edge zero y_ = 1. In the region X4 < x < 
1 — Q, where 

Xi = TTVQ' (75) 

the locus consists of the loop 2/2(0) and the line yi(ip) on the real axis between y* and y_ 
(2/2(0) < y* < y~ < 1). The loop j/2 meets with the line y\ at the point y*, and also cuts the 
real axis at the point 2/2(0) (= 0). At x — x 4 , the line Q\{(f>) shrinks to the point 

»-=»• = »W = Q + (2-Q)V(? - (76) 
For x < X4, the line disappears, and the only locus is the loop 2/2(0) which cuts the real axis 
at two points 2/2(0) and 2/2( 7r ) (0 = 2/2(0) < 2/2( 7r ) < 1)- At Q — 0, all the Fisher zeros lie 
only on the loop 2/2(0) for a; < 1. 



D. One-state (Q=l) Potts model 

At Q = 1, the partition function becomes 

Z N (a,x,Q = l) = (ax) N . (77) 

Therefore, all the Fisher zeros lie on the point a = y~ x = 0. On the other hand, in the limit 
Q — > 1, the loop 2/2(0) is reduced to 

y 2 {<f>) = l-xe- i *, (78) 

and the Fisher zeros are uniformly distributed on this circle. The circle 2/2(0) cuts the real 
axis at two points 1 + x and 1 — x. The point 1 + x may be an antiferromagnetic transition 
point, whereas the point 1—x can correspond to a ferromagnetic transition point for x < 1. 
The locus 2/1 ("0) also becomes the identical circle with the locus 2/2(0)- 
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IV. POTTS ZEROS 



From the equation (8) the locus of the Potts zeros is obtained to be 



gi cos^) + iJg 2 ' 



Qi(^) = l-a + 

where g± — (a — l)x and g 2 = g\ sin 2 ip + x — 1. The edge zeros of Qi(ip) are given by 



(79) 



Q± = (a-2)(x-l)±2<<J(l-a)x(x-i) (80) 
from Q\(n) and Qi(0). If g\ > and g 2 > or g 1 < and g 2 < 0, it is easily verified that 

|Qi(V>)-(l-a)| = \ax-l\. (81) 
From Eq. (81) we see that the zeros of Qiijp) lie on a circle 

QiW=Qc + Se»W (82) 
in the complex Q plane with center 

Q c = l-a (83) 

and radius 

r=\S\, (84) 

where 

S = ax-1. (85) 

The argument 6(ip) is given by 

cos0(VO = 2^— ^cos^-l. (86) 
ax — 1 

The one point of the circle Qi(ip), 

Qi(l> = \ or^ = 7r) =2-a(x + l), (87) 

always lies on the real axis. However, if gi > and g 2 < or g\ < and g 2 > 0, the zeros 
of Qi^) lie on the real axis. 

On the other hand, the equation (12) gives the second locus of the Potts zeros 



W) = (a- l)(» + e «2 -(■» + ■>- 2)<* (88) 
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The two points of ^2(0), 

g 2 (0) = 1 (89) 



and 



Qafr) = - 2a, (90) 



can be lie on the real axis. From |A±| = |Aq| we also obtain 



= (. -!)(!-,) 

X 



A. Three special cases: x = 1, x = 0, and a = 1 

In the absence of an external field (H q = or x = 1), the locus Qi{ip) becomes 

Qi(V>) = 1 - a + (a - l)e* e , (92) 

where the argument 6(ip) is simply given by 

0(V>) = 2^>. (93) 

At the same time, the locus Q2(<f>) reduces to 

Q 2 (0 = (a-l)( e *-l), (94) 

which is the identical circle to the locus Qi(ip). The equation (93) means that the Potts 
zeros are uniformly distributed on this circle. In particular, the Potts zeros lie on the unit 



circle at a = 



691 ]. The circle Qi(ip) or Q2{4>) cuts the real axis at two points Q = 



and 2(1 — a). For ferromagnetic interaction J > the physical interval is 

1 < a < 00 (00 > T > 0), (95) 

whereas for antiferromagnetic interaction J < the physical interval 

< a < 1 (0 < T < 00). (96) 

It should be noted that the point 2(1 — a) lies on the positive real axis for the Potts an- 
tiferromagnet. For | < a < 1 the point 2(1 — a) locates in < 2(1 — a) < 1, while for 
< a < I it does in 1 < 2(1 — a) < 2. The locus of the Potts zeros consists of the circle and 
the isolated point at Q — 1 for < a < ^. 
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In the limit H q — > — oo (x — > 0), the partition function of the Q-state Potts model is 
again given by Eq. (60), and the symmetry of the Hamiltonian is that of the (Q — l)-state 
Potts model in zero external field. At x — 0, the locus Qi^) (Eq- (88)) of the Q-state Potts 
model reduces to 

Q 2 (0) = 2 - a + (a - l)e^, (97) 

which is a circle with center Q c = 2 — a and radius r = |a — 1|. The equation (97) is also 
obtained by replacing Q in Eq. (94) with Q — 1. The circle Q2{4>) cuts the real axis at two 
points Q — 1 and 3 — 2a. For | < a < 1 the point 3 — 2a locates in 1 < 3 — 2a < 2, while for 
< a < ^ it does in 2 < 3 — 2a < 3. The locus of the Potts zeros consists of the circle Q2(4>) 
and the isolated point at Q = 2 for < a < \. The isolated point Q = 2 is also obtained by 
replacing Q in the isolated point Q — 1 for x — 1 with Q — 1. 

At a = 1 (T = oo), the eigenvalue A = a — 1 disappears. Therefore, the Potts zeros are 
completely determined by the locus Qi(ip). All the Potts zeros lie on the point 

Qi(^) = l-x (98) 

for a = 1. 



B. a > 1 and x > 1 

In this case #i and #2 are always positive and the Potts zeros of the ferromagnetic model 
(ferromagnetic Potts zeros), Qi(ip), lie on a circle where the eigenvalues satisfy 

M = TEEjji >t (99) 

|A | V a-1 

which implies that the locus Q2{4 > ) d° es n °t appear. The circle Qi(ip) always cuts the real 
axis at the point Qi(f ) (< 0). Figure 8 shows the loci of the Potts zeros in the complex Q 
plane for a = 3 which are centered at 

Q c = -2. (100) 

The radius of the circle Qi(ip) is 

r = 2 (x = 1), 5 (x = 2), 8 (x = 3) (101) 
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for a = 3. For x — 1 two edge zeros Q± cut the real axis at the origin, whereas they move 
away from the origin as x increases. The Potts edge zeros are determined by the edge angle 

(a-2)x + l 



9 = 6(iP = 0) = cos" 1 
and for a = 3 the edge zeros are located at 



ax — 1 



(102) 



and 



Q± = 0(x=l), 



Q ± = l±Ai (x = 2), 



Q± = 2±4V3i (x = 3). 



On the circle Qi(ip), the density of zeros g(9(tp)) is given by 



sin ■ 



2vr Jsin 2 | - sin 2 f 



(103) 



(104) 



(105) 



(106) 



for \6\ > 9 , and 

9(0) = (107) 

for \9\ < 6q. For x — 1, — 0, and the Potts zeros are distributed on a circle with the 
uniform density of zeros 

9(0) = ^- (108) 
However, for x > 1 the density of zeros at the Potts edge zeros Q± diverges, that is, 

1 



9(0 ~ ) 



y/Q(9) - Q±(0 ) 



(109) 



In this case, the Potts edge zero can be called the Potts edge singularity with the edge 
critical exponent ji e (= — |), and the density of zeros is expressed as 



9(0~9 )~(Q-Q ± )^. 



(110) 



C. a > 1 and x < 1 



For di < a < a 2 , where 



3x + 1 

(x + iy 



(in) 
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and 

a 2 = — !— , (112) 
1 — x 

the locus of the ferromagnetic Potts zeros consists of the loop Q 2 (0) {<P* < 4> < 27r — 0*) and 
the line Qi(tp) (V** < ^ < vr) on the positive real axis between Q* = Qi(ip*) an d Q+ — Qi(Q) 
(0 < < The loop Q 2 meets with the line Q\ at the point 

Q* = Q 2 (0*) =Q 2 (27T-^), (113) 

where |A + | = |A_| = |A |. The loop Q 2 ((f)) cuts the real axis at two points Q 2 (re) and 
(> ^2(71"))- The sign of Q 2 (7r) is positive for a\ < a < 03, where 

" 3 = 2(1+1)' (H4) 

and negative for a 3 < a < a 2 . Figure 9 shows the locus of the Potts zeros for a = § and 
x = I from which we obtain 

Q 2 (7r) = ~, Q* = \, Q + = J + = 0.957. (115) 

At a = ai, 0* = 7T, the other edge zero Q_ appears, and the loop Q 2 {4>) shrinks to the 
point 

Q_ = Q. = Q 2 fr) = (|^|) 2 . (116) 

For a < ai, the only locus is the line Q\(jp) on the positive real axis between Q_ and Q + 
(0 < Q- < Q+). On the other hand, at a = a 2 , 0* = 0, and the line Qi(ip) shrinks to the 
point 

Q+ = Q«. = Q 2 (0) = 1. (117) 

For a > a 2 , all the Potts zeros lie on the loop Q 2 (0) which again cuts the real axis at two 
points Q 2 (n) and Q 2 (0) (Q 2 (7i) < <? 2 (0) = 1). The sign of Q 2 {n) is positive for a 2 < a < 03 
and negative for a > a 3 . 

D. a < 1 and x > 1 

Because <?i is always negative, the zeros of Qi(ip) lie on a circle if g 2 < (-00 < "0 < t 1 " - V'o) 
and on the real axis if g 2 > (0 < -0 < -0 O or t 1 " — V'o < ^ < tO- F° r < a < a 4 , where 

«. = ; i T , (118) 
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the locus of the Potts zeros of the antiferromagnetic model (antiferromagnetic Potts zeros) 
consists of the line Qi(ip) (0 < ip < ipo and ir — -0 O < i> < i>*) on the real axis between 
Q- = Qi(0) (< 0) and Q* = Qi(ip*) (> 0), the circle Qi(ip) (ip < i> < ^ - *0o), and the 
loop Qi{$) (—0* < < 0*), inside the circle Qi- The loop Q 2 meets with the line Q\ at the 
point 

Q* = Q 2 (0*) = Q 2 (-0*), (119) 

where |A + | = |A_| = |A |. The circle Qi(ip) cuts the real axis at two points Q\ (> 0) and 
<2i(f) (1 < <5i(f) < 2), where the points Q\ and Qi(f ) are given by 

Qi = Qi(4>o) = Qi(* - Vo) = Qc + S = a(x - 1) (120) 

and 

Ql GD =Qc-S = 2-a(x + l), (121) 

respectively. Similarly, the loop also cuts the real axis at two points and Q 2 (0) (0 < 
Q* < Q 2 (0) = 1). For example, figure 10 shows the locus for a = and x — 2. In this case 
we obtain 

1Q fi 1 Q /7r\ 17 

Q_ = — - - -p= = -4.583, Q° = — , Q* = — , Qi - = — , (122) 
V 10 ' Vl 10' V 20' V \2J 10' 1 ; 

and the center Q c and the radius r for the circle 

9 4 

Q c = — , r = -. (123) 
^ 10' 5 V ; 

At a = a 4 , two points and Q° 011 the real axis meet at 

Q* = Qi = (124) 
other three points (?i(§), Q 2 (0) and Q 2 (7r) on the real axis also meet at 

Qi(|) = Q 2 (0) = Q 2 (tt) = 1, (125) 
and the circle Qi(ip) and the loop Q 2 (0) become the identical locus as a circle with center 

Qc = -4t ( 126 ) 

rr + 1 

and radius 

r = 7TT' < 127 > 
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On this circle three eigenvalues have the same magnitude 

|A + | = |A_| = | A | = -4t- (128) 
x + 1 

Therefore, the locus consists of the line Qi{ip) on the real axis between Q- and Q* and 
the circle. In the region a 4 < a < a\, the circle Qi(ip) disappears, and the line Qi{ip) 
(0 < ip < ip*) on the real axis between Q_ (< 0) and (> 0) again meets with the loop 
^2(0) (0* < < 27r — 0*) at the point Q*. The loop cuts the real axis at two points and 
Q 2 (tt) (0 < Q* < Q 2 (tt) < 1). 

At a — ai, the other edge zero Q + appears, and the loop (^(tt) shrinks to the point 

Q+ = Q* = Q2M = (^—^) 2 - (129) 



,x + I, 

For a > ai, the loop disappears, and the only locus is the line Qi^) on the real axis between 
Q- (< 0) and Q + (> Q-). The sign of Q + is positive (negative) if a < a 5 (a > a 5 ), where 

In the limit x — > 00, the boundary curves a = ai and a = a 4 approach the line a = 0, and 
the Potts zeros lie on the line Qiijp) which approaches the point —00. 



E. a < 1 and x < 1 



For a < 1 and x < 1, gi and #2 are always negative, and the antiferromagnetic Potts zeros 
of the locus Qi(ip) lie on a circle. At a = a 4 , the loci Qi(tp) and Q 2 (0) become identical to 
be a circle on which three eigenvalues again have the same magnitude 

|A+| = |A_| = |A | = -4t. (131) 
x + 1 

This circle cuts the real axis only at the point 

= Q 2 (0) = Q2M = 1. (132) 

For a > 04, all the Potts zeros lie on the circle Qi(V') (0 < ip < 7r) which cuts the real axis 
at the point <5i(f ) (0 < Qi(f ) < 1)- Figure 11 shows the locus of the Potts zeros for a — | 
and rr = which has the center 

Qc = I (133) 
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the radius 

r = |, (134) 

the edge zeros 

Q± = ^:±4\h = 0.140 ±0.379*, (135) 
50 5 V 5 

and the point on the positive real axis 

On the other hand, for < a < a 4 , the locations of Potts zeros are completely determined 
by the loop Q2{4>) (0 < < 27r) which cuts the real axis at two points ^(0) (= 1) and 
Q2(tt) (1 < Q2(tt) < 3). The locus Q2{4>) nas a crescent-like shape near a = a 4 , but an 
egg-like shape far away from a = a 4 . Figures 12 and 13 show the loci of the Potts zeros for 
a = \ and x — \ and for a = and x = ^, respectively. Figure 12 shows a crescent-like 
shape with 

Q 2 (7r) = |, (137) 
while figure 13 has an egg-like shape and 

144 

Q 2 (vr) = — . (138) 

For a e < a < a 4 , where 

a 6 = } ~ X . , (139) 
6 2(x + l)' v ; 

the point (^(tt) locates in 1 < Q2(n) < 2, while for < a < a6 it does in 2 < Q2(tt) < 3. 

The locus of the Potts zeros consists of the loop Q2{4>) an d the isolated point at Q = 2 for 

< a < a 6 . 



V. CONCLUSION 



We have investigated the properties of the Fisher zeros in the complex temperature plane 
and the Potts zeros in the complex Q plane of the Q-state Potts model in an arbitrary nonzero 
external magnetic field, finding two master equations to determine the partition function 
zeros from the exact partition function of the one-dimensional model. The distribution of 
the Fisher zeros or the Potts zeros is determined by the interplay between these master 
equations. We have discussed the Fisher zeros of the Potts model in an external magnetic 
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field for any real value of Q > 0. We have also studied the Potts zeros in the complex Q 
plane of the ferromagnetic and antiferromagnetic Potts models for nonzero magnetic field. 

Some circle theorems have been found for both the Fisher zeros and the Potts zeros. All 
Fisher zeros lie on a circle for Q > 1 and x > 1 except Q = 2 (Ising model) whose zeros lie on 
the imaginary axis. All Fisher zeros also lie on a circle for any value of Q when x = 1 (except 
Q — 0, 1 and 2) or x = (except Q — 1, 2 and 3). All Potts zeros of the ferromagnetic 
model lie on a circle for x > 1. When x = 1 or x = 0, all Potts zeros lie on a circle for both 
the ferromagnetic and antiferromagnetic models. All Potts zeros of the antiferromagnetic 
model with x < 1 also lie on a circle for (x + l) -1 < a < 1. It has been found that a part of 
the Fisher zeros or the Potts zeros lie on a circle for the specific ranges of x. 

We have shown that some Fisher or Potts zeros can cut the positive real axis. Further- 
more, the Fisher zeros or the Potts zeros lie on the positive real axis for the specific ranges 
of x. We have also calculated and discussed the density of zeros. The density of zeros at the 
Fisher edge singularity diverges, and the edge critical exponents at the singularity satisfy 
the Rushbrooke scaling law. We have found the Potts edge singularity in the complex Q 
plane which is similar to the Fisher edge singularity in the complex temperature plane. 
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Re(y) 

FIG. 1: The locus of the Fisher zeros in the complex y = e _/3J plane of the three-state Potts model 
for x = e l3Hq = 1, 2 and 3. For comparison, the zeros of a finite-size system (N = 100) are also 
shown (open circles for x = 1, pluses for x = 2, and open triangles for x = 3, respectively). 
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FIG. 6: The locus of the Fisher zeros of the Q = \ state Potts model for x = 3. The zeros on the 
real axis lie between — oo and y+ = —0.506 and between = 4 and oo. Most of them are omitted 
in the figure. 
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FIG. 8: The locus of the Potts zeros in the complex Q plane for a = 3. For comparison, the zeros 
of a finite-size system (N = 100) are also shown (open circles for x = 1, pluses for x = 2, and open 
triangles for x = 3, respectively). 
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